Analysis and control of network synchronizability 
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Abstract. In this paper, the investigation is first motivated by showing two examples of 
simple regular symmetrical graphs, which have the same structural parameters, such as average 
distance, degree distribution and node betweenness centrality, but have very different synchro- 
nizabilities. This demonstrates the complexity of the network synchronizability problem. For a 
given network with identical node dynamics, it is further shown that two key factors influencing 
the network synchronizability are the network inner linking matrix and the eigenvalues of the net- 
work topological matrix. Several examples are then provided to show that adding new edges to a 
network can either increase or decrease the network synchronizability. In searching for conditions 
under which the network synchronizability may be increased by adding edges, it is found that for 
networks with disconnected complementary graphs, adding edges never decreases their synchro- 
nizability. This implies that better understanding and careful manipulation of the complementary 
graphs are important and useful for enhancing the network synchronizability. Moreover, it is found 
that an unbounded synchronized region is always easier to analyze than a bounded synchronized 
region. Therefore, to effectively enhance the network synchronizability, a design method is finally 
presented for the inner linking matrix of rank 1 such that the resultant network has an unbounded 
synchronized region, for the case where the synchronous state is an equilibrium point of the network. 

Keywords. Complex network, Complementary graph, Synchronizability, Synchronized region, 
Inner linking matrix. 

1 Introduction and problem formulation 

The subject of network synchronization has recently attracted increasing attention from various 
fields (see [B[2l[3lll[ia[H[Il[22l[23l[2l[25l[28]and references therein) . Of particular importance 
is how the synchronizability depends on various structural parameters of the network, such as aver- 
age distance, clustering coefficient, coupling strength, degree distribution and weight distribution, 
among others. Some important results have been established for such concerned problems based on 
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the notions of master stability function and synchronized region [T] [TT| [l"4l [TBI fl8| 129] . Some inter- 
esting relationships between synchronizability and structural parameters of networks have also been 
reported, e.g., smaller average network distance does not necessarily mean better synchronizability 
|17j . therefore the betweenness centrality was proposed as a good indicator for synchronizability 
[8]. And two networks with the same degree sequence were constructed in a probabilistic sense to 
demonstrate that they can have different synchronizabilities [26], showing that synchronizability 
has no direct relations with degree distribution. Moreover, the effect of perturbations of coupling 
matrices on the synchronizability was studied in [27] . Motivated by all these research works, this 
paper attempts to further explore the analysis and control problems of synchronizability for various 
complex dynamical networks. 

Consider a dynamical network consisting of N coupled identical nodes, with each node being an 
n-dimensional dynamical system, described by 

N 

Xi = f(xi) - c} j a i jH(xj), i = 1,2, ••• ,N, (1) 
j'=i 

where rcj = (xn, Xi2, ••• , £m) € is the state vector of node i, /(•) : M. n — > W 1 is a smooth vector- 
valued function, constant c > represents the coupling strength, H(-) : M n — > W 1 is called the inner 
linking function, and A = (ay)jv><jv is called the outer coupling matrix or topological matrix, which 
represents the coupling configuration of the entire network. This paper only considers the case that 
the network is diffusively connected, i.e., A is irreducible and its entries satisfy 

N 

au = - aij, i = 1,2, ■■■ ,N. 

Further, suppose that, if there is an edge between node i and node j, then = aji = —1, i.e., 
A is a Laplacian matrix. Therefore, is an eigenvalue of A with multiplicity 1, and all the other 
eigenvalues of A are strictly positive, which are denoted by 

= Ai < A 2 < A 3 < • • • < X N . (2) 

The dynamical network (pQ) is said to achieve (asymptotical) synchronization if 

x\{t) — > x%(t) —►•••—> XN(t) —f s(t), as t — ► oo, (3) 

where, because of the diffusive coupling configuration, the synchronous state s(t) € M. n is a solution 
of an individual node, i.e., s(t) = f(s(t)). Here, s(t) can be an equilibrium point, a periodic orbit, 
or even a chaotic orbit. 

As shown in [141 118] . the local stability of the synchronized solution x\(£) = X2(t) = ■■■ = 
%N(t) = s(t) can be determined by analyzing the following so-called master stability equation: 

dj=[Df(s(t)) + aDH(s(t))]uj, (4) 

where q€1, and Df(s(t)) and DH(s(t)) are the Jacobian matrices of functions / and H at s{t), 
respectively. 



The largest Lyapunov exponent L max of network (pQ), which can be calculated from system (|4|) 
and is a function of a, is referred to as the master stability function. In addition, the region S 
of negative real a where L max is also negative is called the synchronized region of network ([1]). 
Based on the results of [14, 18J, the synchronized solution of dynamical network (pQ) is locally 
asymptotically stable if, and only if, 

c\ k £S, A; = 2, 3, ■ ■ ■ ,N. (5) 

The synchronized region S can be an unbounded region, a bounded region, an empty set, or a 
union of several such regions. 

Obviously, for given node dynamics of a linearly coupled network, two key factors influencing 
the synchronizability are the inner linking matrix H (•) = H and the eigenvalues of the topological 
matrix A. The inner linking matrix is directly related to the synchronized region, as studied 
in [Tl~| [T2l [TBI I20j - The larger the synchronized region, the easier the synchronization. The 
topological matrix, on the other hand, is directly related to condition ©. If S is an unbounded 
sector (— oo, a], the eigenvalue A2 of A determines the synchronizability [23]; if S is a bounded 
sector [01,02]) the ratio r(A) = determines the synchronizability [JJ. No matter what the 
synchronized region is, the larger the A2 and r(A) are, the easier the synchronization is. This paper 
will further study this issue with more careful analysis. 

The rest of this paper is organized as follows. In Section 2, two simple graphs on six nodes are 
given to show that networks with the same structural parameters, such as average distance, degree 
distribution and betweenness centrality, can have different synchronizabilities. If the synchronized 
region S is unbounded, adding edges never decreases the synchronizability, but this may not be 
true if S is bounded. In Section 3, a class of networks with disconnected complementary graphs 
are discussed. For such networks, adding edges never decreases the synchronizability no matter 
what type of region S is. In Section 4, an inner linking matrix of rank 1 is designed for realizing 
unbounded synchronized regions in the case that the synchronous state is an equilibrium point. In 
Section 5, some network synchronization examples are provided to illustrate the theoretical results. 
The paper is concluded by the last section. 

Throughout this paper, for any given undirected graph G, eigenvalues of G mean eigenvalues 
of its corresponding Laplacian matrix. Notations for graphs and their corresponding Laplacian 
matrices are not differentiated, and networks and their corresponding graphs are not distinguished, 
unless otherwise indicated. 

2 Two simple graphs that tell the main idea 

In [26], two networks with the same degree sequence were constructed in a probabilistic sense to 
show that they can have different synchronizabilities. In this section, the two simple graphs G\ 
and G2 on six nodes, shown in Figs. 1 and 2, are considered, where G\ is a typical bipartite graph 
with many interesting properties. 
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Fig. 1 Graph G x Fig. 2 Graph G 2 

Obviously, graphs G\ and G 2 have the same degree sequence, where the degree of every node is 
3; the same average distance |; and the same node betweenness centrality 2 [TIE]. Although these 
two graphs have the same structural characteristics, their corresponding networks have different 
synchronizabilities, as shown below. Their Laplacian matrices are 

1 1 A 


1 

1 ' 
-3 1 
1 -3 / 

respectively. The eigenvalues of G\ are 0, 3, 3, 3, 3 and 6; the eigenvalues of G 2 are 0, 2, 3, 3, 5 
and 5. Obviously, A 2 (G X ) = 3 > A 2 (G 2 ) = 2, and r(Gi) = 0.5 > r(G 2 ) = 0.4. Therefore, the 
synchronizability of network G\ is better than that of network G 2 . 

Graphs G\ and G 2 have the same structural parameters, but it is clear that they have different 
average clustering coefficients, denoted by C(Gi),i = 1,2, with G(G 2 ) > C{G\). As mentioned 
above, the clustering coefficient does not have direct relation to synchronization [T7]. For example, 
globally coupled graphs have the largest clustering coefficient, 1, and they have the best synchro- 
nizability. However, for the above two graphs, the larger clustering coefficient does not indicate 
better synchronizability. This is demonstrated by the following process of adding edges. 

Consider enhancing A 2 and r by adding edges to G 2 . For this purpose, the following lemma is 
needed. 

Lemma 1 [15] For any given connected undirected graph G of size N, its nonzero eigenvalues 
indexed as in ([2]) grow monotonically with the number of added edges, that is, for any added edge 
e, Xi(G + e) > Xi(G), i = 1, ■ ■ • ,N. 

By Lemma 1, obviously if the synchronized region is unbounded, adding edges never decreases 
the synchronizability. However, for bounded synchronized regions, this is not necessarily true. 
For example, adding an edge between node 1 and node 3 in graph G 2 (Fig. 2), denoted by 
e{l, 3}, leads to a new graph G 2 + e{l, 3}, whose eigenvalues are 0, 2.2679, 3, 4, 5 and 5.7321. Thus, 
r(G 2 + e{l,3}) = 0.3956 is even smaller than r(G 2 ) = 0.4. This means that the synchronizability 
of network G 2 + e{l,3} is worse than that of network G 2 . Adding a new edge between node 1 
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and node 4 instead, one gets r(G 2 + e{l,3}) < r(G 2 + e{l,3} + e{l,4}) = 0.3970 < r(G 2 ). This 
means that, the synchronizability of network G 2 + e{l,3} + e{l,4} is better than G 2 + e{l,3}, 
but still worse than G 2 . Therefore, by adding edges, the network synchronizability may increase 
or decrease, for which no general rule has been found to date. 

On the other hand, during the process of adding edges, average distance decreases and average 
clustering coefficient increases. But this does not indicate better synchronizability, consistent with 
the conclusion in |17j . 

It was shown [8j that the synchronizability is always improved as the maximum betweenness 
centrality is reduced, which is consistent with the conclusion of [T7]. In the above two graphs, 
however, it shows that the same betweenness centrality does not necessarily mean the same syn- 
chronizability. On the other hand, adding three edges between nodes 1 and 6, 2 and 3, 3 and 4, 
respectively, in graph G±, and then computing their corresponding eigenvalues, it can be verified 
that the networks corresponding to the resulting graphs G = G\ + e{l, 6} + e{2, 3} + e{3, 4} and 
Gi have the same synchronizability. However, in this case, the maximum betweenness centrality of 
G, is smaller than that of Gi, 2. This shows that the smaller betweenness centrality does not 
necessarily indicate better synchronizability, revealing the complexity in the relationship between 
synchronizability and network structural parameters. 

Note that adding edges in G\ also increases the clustering coefficient and decreases the average 
distance, but this does not result in the increase of synchronizability. In the following, it explains 
why adding three edges in G\ does not increase the synchronizability. 

3 Networks with disconnected complementary graphs 

For a given graph G, the complement of G is the graph containing all the nodes of G, and all the 
edges that are not in G. The complementary graph of G is denoted by G c . For example, the com- 
plementary graphs of G% and G 2 in Fig. 1 and Fig. 2 are shown in Fig. 3 and Fig. 4, respectively. 
In the previous section, it shows that adding edges sometimes decreases the synchronizability. How- 
ever, for a class of graphs with disconnected complementary graphs, this never occurs. In order to 
discuss such networks, the following lemma is needed. 
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Fig. 3 Graph G\ Fig. 4 Graph G% 

Lemma 2 [15] For any given graph G, the following statements hold: 



(i) Ajv(G), the largest eigenvalue of G, satisfies \n(G) < N. 

(ii) Xn(G) = N if, and only if, G c is disconnected. 

(iii) If G c is disconnected and has (exactly) q connected components, then the multiplicity of 
X N (G) =N is q- 1. 

(iv) Xi(G c ) = N - Ajv_ i+2 (G), 2<i<A. 

The complementary graph of G\ (Fig.l) is shown in Fig. 3, which is disconnected. The largest 
eigenvalue of G\ is 6, which remains the same when the graph receives additional edges. Hence, 
combining with Lemma 1, the synchronizability of the networks built on graph G\ never decrease 
with adding edges. Although this is true, adding any 3 edges to graph G\ does not enhance the 
synchronizability, since the least nonzero eigenvalue A 2 = 3 of G\ has multiplicity 4 (the multiplicity 
of the largest eigenvalue in G\). This is due to the fact that, for any graph G, rank(Aj/ — (G + e)) < 
rank(A,/ - G) + 1. 

According to Lemma 2, the multiplicity of the largest eigenvalue of a graph G is related to 
the number of connected components of its complement G c . Hence, with the same number of 
edges, generally the synchronizability of the networks built on graph G is better when G c has two 
connected components than the case that G c has more than two connected components. This is 
because in the latter case, the multiplicity of the largest eigenvalue of G is larger than 1, i.e., 
some edges have no contributions to the synchronizability. On the other hand, the multiplicity of 
the largest eigenvalue of G c (i.e., the multiplicity of the least nonzero eigenvalue of G) should be 
large in order to reduce the number of edges needed to enhance the synchronizability. Therefore, 
better understanding and careful manipulation of complementary graphs are useful for enhancing 
the network synchronizability. And, at least for dense networks, the complementary graphs are 
easier to analyze than the original graphs, e.g., G\ in Fig. 3 is simpler than G\ in Fig. 1. 

Remark 1 The graphs shown in Figs. 1 and 2 can be generalized to graphs of size N = In. 
Suppose that graph G\ is bipartite in the sense that it contains two sets of nodes, each set containing 
n isolated nodes, and each node in one set connects to all the nodes in the other set, i.e., the 
complementary graph of G\ is two separated fully connected subgraphs of size n. Graph G 2 is 
composed of two fully connected subgraphs and n edges connecting each node in one subgraph 
to the corresponding node in the other subgraph. In this case, the least nonzero and maximum 
eigenvalues of G\ are y and N, respectively, and r(G\) = |. On the other hand, the least nonzero 
and maximum eigenvalues of G 2 are 2 and y + 2 [15], respectively, with r(G 2 ) = — > 

as N — > +oo. Therefore, these two graphs have the same structural parameters but have very 
different synchronizabilities. 

4 Designing the inner linking matrix 

From the above two sections, it can be seen that the bounded synchronized regions are more 
complicated than the unbounded synchronized regions. Thus, the synchronizability is easier to 
analyze when the synchronized region is unbounded. Hence, it is interesting to find out how to 
design the inner linking matrix such that the network synchronized region is unbounded. 



If the synchronous state is an equilibrium point, then both Df(s(t)) and DH(s(t)) in (|4]) reduce 
to constant matrices, denoted by F and H , respectively. In this case, system (j3J) becomes 

oj = [F + aH]u. (6) 

Hence, the synchronized region S becomes the stability region of F + aH with respect to parameter 
a. In this section, consider the design of an H such that F + aH has an unbounded stable region. 
It is well known that if H is an anti-stable matrix (e.g., H = I n ), F + aH has an unbounded stable 
region [12] . However, if H is of full rank, it means that the coupling in the network is a full state 
coupling among nodes, so the coupling cost may be high. For this reason, consider the design of 
an H of rank 1 such that the stability region for F + aH is unbounded. In this case, the coupling 
can be viewed as an input-output coupling as in control systems [S], or an observer-based coupling 

ma. 

Theorem 1 Given a matrix F G R raxn , there exists a matrix H € R nxn of rank 1 such that 
the stability region of F + aH with respect to parameter a contains (— oo, ai], a\ < 0, if and only 
if every unstable eigenvalue of F is corresponding to only one Jordan block. 

Proof Without loss of generality, suppose a\ = —1. If F + aH is stable and H is of rank 1, it 
means that (F, H) is stabilizable, so that every unstable eigenvalue of F must be corresponding to 
only one Jordan block. 

On the other hand, if every unstable eigenvalue of F is corresponding to only one Jordan block, 
one may take a column vector b such that (F, b) is stabilizable. Then, there exists a row vector k 
such that F — bk is stable, i.e., there exists a matrix P = P T > such that 

(F — bk)P + P(F — bk)' < 0. 

Let kP = y, so that the above inequality becomes 

FP + PF' -by- y'b' < 0. (7) 

By the canonical projection lemma in H^ control theory [9], there exists y such that the above 
inequality holds if, and only if, there exists a scalar f3 > such that 

FP + PF' - (3bb' < 0. (8) 

Since P is a matrix to be determined, in the above inequality, suppose j3 = 2 without loss of 
generality. Then 

FP + PF' - 2bb' < 0. (9) 

Obviously, if (F, b) is stabilizable, then there exists P = P T > such that ([!]) holds. And, when 
holds, for any > 2, (jHJ) holds. Take y = b' , i.e., k = b'P^ 1 . By the above inequalities, one 
knows that F — abk is stable for all a E (-co, —1]. Therefore, H = bk is the matrix to be found. □ 
Remark 2 Let Z\ = kxi and the inner linking function H(xj) = bkxj in network ([T]). Then Z{ 
can be viewed as the output of node i of (pQ) and the linking function bkxj can be viewed as the 
influence of the output of node j to the other nodes. Clearly, the above coupling is simpler than 
full state couplings. 



Remark 3 If F is stable, i.e., the node system is locally stable, then there is always a matrix 
H of rank 1 such that the resulting network has an unbounded synchronized region, as shown by 
the examples given below. In addition, if there exists an unstable eigenvalue of F whose number of 
Jordan blocks is more than 1, one may take an inner coupling matrix H such that its rank equals 
the maximum number of Jordan blocks among the unstable eigenvalues of F. On the other hand, 
one may design an H such that F + aH has an unbounded stable region. Moreover, one may also 
design an H such that the unstable region of F + aH is unbounded, if desired, which is useful for 
desynchronization problems. 



5 Examples 



Example 1 Consider the network (fT]) consisting of the third-order smooth Chua's circuits [2T] , in 
which each node is described by 



An = —kaxn + kaxi2 — ka(axf 1 + bxn), 
x i2 = kxn - kx i2 + kx i3 , 
Xi3 = -k/3x i2 - kjx i3 . 

The vector x, in (fT]) is (xn, Xi 2 , Xi%) T here. Linearizing (|1U|) at its zero equilibrium gives 



(10) 
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(11) 



Take k = l,a = — 0.1,/? = —1,7 = l,a = 1,6 = —25. Then F is stable, i.e., the node system 
10]) is locally stable about zero. Further, take the inner linking matrix [5] 

( 0.8348 9.6619 2.6591 \ 
0.1002 0.0694 0.1005 
-0.3254 -8.5837 -0.9042 



H 



\ 



) 



Then, by simple computation, one knows that F + aH has two disconnected stable regions: S\ = 
[—0.0099, 0] and £2 = [—2.225, —1). Therefore, the entire synchronized region is U ^ • Moreover, 
suppose that the number of nodes is N = 6, and the outer coupling matrix A is equal to the G\ 
in Section 2. According to the eigenvalues of Gi given in Section 2, one may take the coupling 
strength c = 7^9. Then, for every eigenvalue of G\, one has cAj € S 2 - By ([5]), network ([T]) specified 
with the above data achieves synchronization. However, for the outer coupling matrix G 2 given in 
Section 2, for any coupling strength c € [0.002, +00), © does not hold. Therefore, for the above 
node equation, inner coupling matrix and coupling strength, the network built on graph G\ in Fig. 
1 achieves synchronization, but the network built on G 2 in Fig. 2 does not synchronize. Figs. 5 
and 6 show their synchronous and non-synchronous behaviors, respectively. 



,0 20 30 40 50 60 70 80 _0 10 20 30 40 50 60 70 80 

Fig. 5 Network on graph G\. Fig. 6 Network on graph G*2. 

From Figs. 5 and 6, one can see that different topological structures result in very different 
dynamical behaviors for the smooth Chua's circuit network. 

Example 2 Consider designing the inner coupling matrix such that the corresponding network 
has an unbounded synchronized region, where the node system is as in Example 1. In this case, 
F is stable, so for any column vector b, (F, b) is stabilizable, e.g., b = (0, 0, 1) T . By the method of 
Theorem 1, one gets k = (0.0708, —0.15590,0.4296). Then, change the inner coupling matrix H in 
Example 1 to H = bk, while keeping the other parameters unchanged. This F + aH is stable for 
a E (— oo, —1]. Consequently, the corresponding network has an unbounded synchronized region. 
Figs. 7 and 8 show their similar synchronization behaviors of the two networks built on graphs G\ 
and C?2, respectively. 




Fig. 7 Network on graph G\. Fig. 8 Network on graph G*2. 

One can see that the network solutions synchronize a little faster in Fig. 7 than that in Fig. 8. 
This also demonstrates that the synchronizability of the network built on Gi is better than the one 
on G2, which verifies the theoretical analysis given in the previous sections. 



6 Conclusion 



In this paper, the synchronizability of complex dynamical networks, which is directly related to the 
inner linking matrix and the topological matrix, has been carefully discussed. Two simple graphs 



have been given to show that networks can have different synchronizabilities even when they have 
the same average distance, node betweenness centrality and degree distribution. It has also been 
shown that the larger clustering coefficient, smaller betweenness centrality and average distance do 
not necessarily imply better synchronizability. This demonstrates the complexity in the relationship 
between the synchronizability and network structural parameters. The most significant discovery of 
this paper is that if the synchronized region is bounded, adding edges can either increase or decrease 
the network synchronizability; however, for networks with disconnected complementary graphs, 
adding edges never decreases their synchronizability. Therefore, better understanding and careful 
manipulation of complementary graphs is useful for enhancing network synchronizability. Moreover, 
unbounded synchronized regions are easier to analyze than the bounded ones. To effectively enhance 
the network synchronizability, a design method for the inner linking matrix of rank 1 is finally 
provided such that the resulting network has an unbounded synchronized region for the case where 
the synchronous state is an equilibrium point of the network. 
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